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Abstract 

Lower bounds for the factors entering the standard notions of shear and 
torsion stiffness for a linearly elastic rod are established in a new and simple 
way. The proofs are based on the following criterion to identify the stiffness 
parameters entering rod theory: the rod's stored-energy density per unit 
length expressed in terms of force and moment resultants should equal the 
stored-energy density per unit length expressed in terms of stress components 
of a Saint-Venant cylinder subject to either flexure or torsion, according 
to the case. It is shown that the shear factor is always greater than one, 
whatever the cross section, a fact that is customarily stated without proof in 
textbooks of structure mechanics; and that the torsion factor is also greater 
than one, except when the cross section is a circle or a circular annulus, a 
fact that is usually proved making use of Saint- Venant's solution in terms 
of displacement components. 

Keywords: Rod theory, shear stiffness, torsion stiffness, shear factor, tor- 
sion factor 



1 Introduction 

When a direct approach is chosen to expound the standard theory of linearly 
elastic rods with a straight axis C, the Principle of Virtual Working is laid 
down: 

/ (T7 + Ne + Mip + M t if) t )dx 3 =: W i = W e := / (jpv + qw + ap + c t ^)dx z . 
J c J c 

(!) 

The internal working W" collects four pairs (T, 7), . . . , (Mt,tp t ) of kinematic 
and dynamic fields in duality, with T, N, M, and Mt the shear force, normal 
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force, bending moment and twisting moment and with 7, e, tp and tp t the 
shear, extension, flexion and torsion measures. The geometric compatibility 
conditions are: 

7 = v + tp, e = w' , tp = tp and tp t = 

where v and w are the transverse and axial displacements, (p and 1? are 
the rotations about, respectively, the a^-axis and the X3-axis (a superscript 
prime denotes differentiation with respect to £3; see Fig. [1]); the external 




Figure 1: 

working W e sets these kinematic variables in duality with the applied loads 
per unit length p,q,c and ctH The constitutive equations are: 

T = s s j,...,Mt =S t ip t , 

all the stiffness moduli s s , . . . ,St being assumed positive; the stored-energy 
density per unit length is: 

1 T 2 1 Ml 

(see, e.g., [I]). 

2 To keep things simple, in laying down flj we have ignored the contributions to external 
and internal workings of concentrated loads, axial pin-junctions et similia. 
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To make use of such a theory — in fact, to make use of any structure 
theory formulated by way of a direct approach — the issue of parameter 
identification must be dealt with. In particular, one would like to gather 
from some suitable equilibrium theory for a rod-like three-dimensional body 
subject to flexure and torsion enough information to choose the 'right' shear 
and torsion stiffness parameters 5 s ,St, that is to say, to choose those param- 
eters so as to guarantee that the qualitative and quantitative predictions of 
the ensuing rod theory allow for a quick, but technically sufficient, evaluation 
of the behavior of the corresponding three-dimensional body. 

Customarily, the three-dimensional theory one picks is the treatment 
of the cases of flexure and torsion of the Saint-Venant Problem in classic 
linearly isotropic elasticity. With this choice, one is led to set: 

S s — , St — , \l) 

Xs Xt 

where G is the shear modulus of the material Saint- Venant's prismatic cylin- 
der is comprised of, A and J Q are the area and the polar moment of inertia 
of the cylinder's cross section, and the dimensionless coefficients Xs an d Xt 
are the shear and torsion factors. Generally, these identifications are justi- 
fied by inspection of the Saint-Venant solutions for flexure and torsion in 
terms of displacement components. Each of these solutions is constructed 
by means of Saint -Venant's semi-inverse method, starting from an educated 
guess about a representation for the displacement field depending on an 
as-short-as-possible list of parameter functions. Since the analytic forms of 
these functions depend solely on the shape of the cross section, the same is 
automatically true for the shear and torsion factors. 

In textbooks of structural mechanics, the shear factor is interpreted as 
a measure of the non-uniformity in the distribution over the cross section 
of Saint-Venant's cylinder of the tangential stresses induced by the shear 
force; it is stated, but rule not proved, that 

X s > 1- (3) 
As to the torsion factor, it is generally stated, and proved, that 

Xt > 1, (4) 

with equality holding if and only if the cross section is a circle or a circular 
annulus. While this facts suggests that (xt ~ 1) measures the deviation from 
central symmetry of the cross section, a precise interpretation is generally 
non attempted; at best, Saint-Venant's conjecture that a circular section 
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minimizes the twist angle among all simply-connected sections with the same 
area is mentioned. 

In this note we approach the issue of parameter identification as is done 
in [2]: we require that a linearly elastic rod stores the same energy per unit 
length as a Saint- Venant cylinder subject to end loads having the same force 
and moment resultant. The expression we use for the rod energy is in terms 
of either the shear resultant force or the twisting moment, that for Saint- 
Venant's cylinder is in terms of stress components. Both to prove (pQ) and 
to characterize the torsion optimality of center-symmetric cross sections, we 
make no use of Saint- Venant's solutions in terms of displacement; the only 
assumptions we borrow from Saint- Venant's problem are those restricting 
the admissible external forces and the class of stress fields. 



2 The shear and torsion factors 

Given a cartesian reference frame {O ; xi, X2> £3} an d the associated or- 
thonormal basis {ei, 62,63}, we consider a Saint-Venant's cylinder whose 
axis is parallel to e% and whose cross section has its barycenter O on the 
X3-axis. The vector x = x\e\ + xi^-i yields the position of a point of a 
typical cross section A with respect to O. On denoting by Sij = 1, 2, 3) 
the cartesian components of the stress tensor, with Saint- Venant we restrict 
attention to cases when both external distance forces and external contact 
forces over the cylinder's mantle are null and, moreover, 

S a ^ = (a, = 1,2). (5) 

Consequently, in order to satisfy the equilibrium equations, not only the 
fields S^a (a = 1, 2) must be independent of X3 but also 

s ■ n = Ssai^a = , x G dA , (6) 

where s = S 3a e a is the cross-sectional traction, n = n a e a is the normal to 
the mantle, and dA is the boundary curve of A. 

When expressed in terms of stress components, the stored-energy density 
for unit volume of an isotropic material is: 

1 1. -a x 



where A and [1 = G are the Lame coefficients, \S\ = (SySy) 1 ^ , and tr S = 
5*11 + S22 + ^33! under Saint-Venant's assumption (J2j) , ([2]) becomes: 

™sv = ^(S 31 + S 32 ) + -S m , E.= — , 
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where E is the Young modulus. Thus, in Saint-Venant's cylinder, the stored- 
energy density splits additively into two terms: the first one is non-null if 
and only if the applied loads induce cross-sectional stress, while the second 
one is non-null if and only if the applied loads induce axial stress. On the 
other hand, the stored-energy density per unit length of a linearly elastic 
rod is: 

III 

2 s7 



when the rod is subjected to a shear force of magnitude T, and is: 

1 Mi 

when the rod is subjected to a torsion moment of magnitude Mt. 

Following [2] , we identify the shear-stiffness parameter by imposing that 

with the integrand on the right side proportional to T 2 ; hence, 



J w sv dA^j = ± j (£% x + Sj 2 ) dA , 



G 



T 1 



USl + Sl 2 )dA 



(8) 



Next, we compare this relation with the first of ([T]): since the factor multi- 
plying the shear modulus G has the dimension of an area, we find it natural 
to give s s the form (JT]) by defining 

Aj A (Sl + Si 2 )dA 



X S ■■= J ^ \ 2 , A = J A dA - ( 9 ) 

Quite similarly, we identify the torsion-stiffness parameter by imposing 



that 



I Mi 



1 



2lT~ 2GJ A {S31 + S32) 



dA. 



with the integrand on the right side proportional to M 4 2 ; hence, 



St 



G 



Mf 



(10) 



This time, the factor multiplying G has the dimension of a moment of inertia; 
to recover the second of (pQ), we set: 

Jof A (Sl + Sl 2 ) dA 



Xt :z 



Jo 



L 



\x¥dA. 



(11) 
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3 Lower bounds 

In this section we prove the lower bounds (pQ) and ([1]) for, respectively, x s 
and Xt- 

3.1 Xs >l. 

Consider the following chain of inequalities: 

A Jj^Sl + S 2 32 ) ^A>aJ S 2 32 dA ^{J A S ^ dA ) 2 ■ ( 12 ) 

Note that the first inequality in (|3.ip holds true with the equality sign if and 
only if S31 is identically null on A; and that the second, which is established 
by making use of Jensen's inequality for convex functions, reduces to an 
equality if and only if S32 is identically constant on 

aE Next, without any 
loss of generality, take the shear force parallel to e 2 , so that 

T= f S 32 dA, 
J A 

(J3JJJ yields: 

A [ (S 2 31 + Si 2 )dA>T 2 , 
J A 

and ([2]) reduces to 

Xs > 1, 

with 

X s = 1 S31 = and S32 = const . 

This set of conditions on the stress field is incompatible with the boundary 
condition ([2]) no matter the shape of the cross section, because the stress 
field must be continuous up to the boundary of the cross section itself. We 
conclude that the strict inequality (P) must hold. 

3-2 xt>l- 
By definition, 

M t = f ( Xl S 32 - x 2 S 31 ) dA. 
J A 

3 A use of this second inequality in similar circumstances is found on p. 475 of [3]- 
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With an application of Pubini's theorem, the denominator in ([2]) can be 
written as 

M ' = L A ( XlS:i2( - x) ~ X2S31 ^) (fi-Mv) ~ toSsiivj) dAdA, 
or rather, equivalently, given that 

/ x iy2 S 32 (x)S 31 (y)dAdA = / Vl x 2 S 32 (y)S 3 i(x) dAdA , 
JAxA J Ay. A 

as 

Ml = [ (x 1 y 1 S 32 (x)S 32 (y)+x 2 y 2 S 31 {x)S 31 (y)-2x 1 y 2 S 32 (x)S 31 (y)) dAdA. 
JAxA v 7 

(13) 

An iterated use of the inequality: 

±2ab <a 2 + b 2 (14) 

yields: 

x m S 32 (x)S 32 (y) dAdA < \ ! (x 2 S 2 32 (y) + ylS 2 32 (x)) dAdA = 

A 1 JAxA v 7 



= / (x 2 S 2 2 (y))dAdA, 

JAxA v 7 

x 2 y 2 S 31 (x)S 31 (y)dAdA <\l (x 2 2 S 2 31 (y) + y 2 S 2 31 (x)) dAdA 

1 JAxA v 7 



AxA 



= I (x 2 Sl(y))dAdA, 
JAxA v 7 

-2 Xl S 32 (x)y 2 S 31 (y)dAdA < [ (x 2 S 2 31 (y) + y 2 S 2 31 (x)) dAdA 

JAxA v 7 

= / (xlSl^ + xlsKy^dAdA. 
JAxA v 7 



AxA 

With this and another use of Fubini's Theorem, it follows from (|3.2|) that 

M 2 < f (x 2 S 2 2 (y) + x 2 S 2 31 (y) + x 2 S 2 31 (y) + x 2 S 2 32 (y)) dAdA = 
JAxA v 7 

= / \\x\\ 2 {Sl 1 (y) + Sl 2 (y))dAdA= ( f \\x\\ 2 dA) [ (S 2 31 (x) + S 2 32 (x)) dA 
JAxA K JA 7 J A 

= Jo f (S 2 1 + S 2 2 )dA, 

J A 



which is tantamount to having from ([2]) that Xt — 1 whatever the shape of 
the cross section^ 

Now, given that equality holds in (|3.2p if and only if a = =F&, Xt = 1 if 
and only if all of the following conditions are identically satisfied in A x A: 

xiS 3 2(y) = ViSvzix), x 2 S 31 (y) = y 2 S 3 i(x), x 1 S 3 i(y) = -y 2 S 32 (x); 

the first two are implied by the last, which can be written as: 

5 3 2(g) _ S 31 (y) _ c 
x% y 2 

or rather, equivalently, 

s = c(-x 2 e 1 +xie 2 ) = ce 3 x x, (15) 

with c a constant. Hence, the cross-sectional traction s must be orthogonal 
to the position vector x all over A up to the boundary, where on the other 
hand it has to be orthogonal to the normal n to satisfy ([2]) . Because of the 
assumed continuity of the stress field up to the boundary, x must then be 
parallel to n all over dA, which is possible if and only if A is a circle or a 
circular annulus. Thus, as anticipated in the introduction, (Xt~ 1) quantifies 
the reduction in torsional stiffness due to the defect in polar symmetry of 
the cross section and the accompanying deviation from the form (13.2|) of the 
cross-sectional traction. 



4 Final remarks 



Remark 1 One can define the extension and bending stiffnesses s e and 5& 
and the relative factors x e an d Xb °f a roc i by posing, respectively, 



for 



and 



liV 2 
2"s7 


2E . 


/ ^33 
'A 


dA 




Aj A Sj,dA 
N 2 






N = 


[ S 33 dA, 

J A 




1 M 2 
2^T 


1 

'~ 2E 


f sh 

J A 


dA 


, Xb 


Jf A S 2 3 dA 
M 2 



(16) 



(17) 



4 A reviewer of our paper remarked that the inequality M t 2 < Jo /^(Sfi + S^) dA can be 
proved by a straightforward application of the Cauchy-Bunyakovsky-Schwarz inequality. 
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for 



M = 



X2S33 dA, J 




A 



J A 



To prove that 



X e > 1 and Xb > 1 



(18) 



is left as an exercise for the interested reader. We note that equality is 
realized in relations (j3J) whatever the shape of the cross section A if and 
only if, respectively, the field S33 is const ant- valued or linear in X2 over A, 
just as it happens to be in Saint- Venant's cases of normal force and pure 
bending. 

Remark 2 The recipe for parameter identification embodied in our defini- 
tions ([2|), ([2]), and for rod stiffnesses is applicable as such to any 
one-dimensional counterpart of a three-dimensional rod-like material body, 
no matter its axis were straight, its cross section constant or its mechan- 
ical response spatially uniform, as is the case for Saint- Venant's prismatic 
cylinder; while linearity in the elastic response is crucial, isotropy is not: all 
those definitions make sense also for transverse isotropy with respect to the 
axial direction. 
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